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ABSTRACT 



Aims. We investigate the nature of 24/jm sources in M33 which have weak or no associated Ha emission. Both bright evolved stars 
and embedded star forming regions are visible as compact infrared sources in the 8 and 24yum Spitzer maps of M33 and contribute to 
the more diffuse and faint emission in these bands. Can we distinguish the two populations ? 

Methods. We carry out deep CO J=2-l and J= 1-0 line searches at the location of 18 compact mid-IR sources and 2 optically selected 
ones to unveil an ongoing star formation process throughout the disk of M33. We use different assumptions to estimate cloud masses 
from pointed observations. We also analyze if the spectral energy distribution and mid-IR colours can be used to discriminate between 
evolved stars and star forming regions. 

Results. Molecular emission has been detected at the location of 17 sources at the level of 0.3 K km s _I or higher in at least one of 
the CO rotational lines. Even though there are no giant molecular clouds beyond 4 kpc in M33, our deep observations have revealed 
that clouds of smaller mass are very common. Estimated molecular cloud masses range between 10 4 and 10 s M , assuming likely 
values of the CO-to-H 2 conversion factor and virial equilibrium. Sources which are known to be evolved variable stars show weaker 
or undetectable CO lines. Evolved stars occupy a well defined region of the IRAC color-color diagrams. Star forming regions are 
scattered throughout a larger area even though the bulk of the distribution has different IRAC colors than evolved variable stars. We 
estimate that about half of the 24 jim sources without an Ho- counterpart are genuine embedded star forming regions. Sources with 
faint but compact Ha emission have an incomplete Initial Mass Function (IMF) at the high-mass end and are compatible with a 
population of young clusters with a stochastically sampled, universal IMF. 
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1. Introduction 

Our knowledge of molecular clouds and of the processes in the 
interstellar medium (hereafter ISM) that lead to the birth of stars 
is mostly based on Galactic studies. Local Group galaxies are 
however sufficiently close to allow individual massive stars and 
molecular clouds to be detected. M33, at a known distance, has 
a high star formation rate per unit area and a low overall ex- 
tinction compared to M31, owing to the moderate gas content 
and low inclination. It is therefore an ideal laboratory for the in- 
vestigation of the relationship of molecular clouds to other ISM 
components and evolutionary scenarios involving blue, low- 
luminosity galaxies. Recent high-resolution optical (HST), in- 
frared (Spitzer; hereafter, IR) and 21 -cm observations (VLA) 
have traced star formation and the ISM throughout the M33 disk 
with high accuracy . Our i nvestigation of the IR emission in M33 
(IVerlev et alJl2007L 120091 1201 d ICorbelli eT^UMk via Spitzer 
high-resolution images has unveiled a variety of star formation 
sites through infrared colors and optical-to-IR ratios. In particu- 
lar, our analysis has shown the existence of two types of IR se- 
lected sources: sources with only diffuse or very faint Ha emis- 
sion, and sources with a definite Ha counterpart. In the former 
sample we can find sites at an early stage of massive star forma- 
tion which give us the opportunity to study individual embedded 
newly born HII regions in a galaxy different than our own. 

Young stellar clusters prior to the phase of gas removal, due 
to photoionization or mechanical force by stellar winds, are em- 
bedded into molecular gas and detectable only at infrared and 



radio wavelengths. Previous searches in M33 have not been suc- 
cessful in detecting embedded clusters. A radio-selected s ample 
of sources in M33 has been analyzed by Buckal ew et al.l {2006) 
who found optically visible counterparts w ith ages between 2- 
10 Myr. Similar ages have been derived by iGrossi et al.l (1201 Oh 
by analyzing the spectral energy distribution of a sample of com- 
pact HII regions. These results point out the paucity of embed- 
ded clusters which might be a short-lived phase of the cluster 
lifetime. A dust abundance lower than usual or a mass spectrum 
of molecular clouds steeper than in our Ga laxy might be respon- 
sible for this result. Buc kalew et al.l d2006l) suggest to analyze an 
IR selected sample, an approach that is now possibl e thanks to 
M33 S pitzer images and 24 /urn source catalogue of IVerlev et al.l 
(120071) . Mid-IR sources which have no visible counterpart in 
the Ho- emission map are generally faint and their nature is 
not obvious. Candidates include evolved clusters, evolved stars 
with dusty envelopes (such as pulsating asymptotic giant branch, 
hereafter AGBs, carbon stars etc.), embedded star forming sites, 
small young clusters which lack massive stars and hence ion- 
izing photons. From the available catalogues we can exclude 
evolv ed clusters (as well as Planetary Nebulae: see IVerlev et al.l 
120071) . since these have already removed their dusty envelope. 
Bright evolved stars have also been catalogued and in this paper 
we discuss the likely presence of contamination. The presence 
of molecular clouds around these sources would instead confirm 
ongoing star formation. 
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A full imaging of molecular clouds compl exes in M33 has 
been completed by the BIMA interferometer (lEngargiola et al.1 
2003b and the FCRAO-14m telescope ( Corbel lil 120031; 



Hev eretal.ll2004 using the 12 CO J=l-0 line. In our Galaxy, 



Giant Molecular Clouds (hereafter GMCs) break up into 
smaller subunits when observed at high spatial resolution 
dRosolowskv et al.l 120031) . However, in M33 both surveys 
(FCRAO and BIMA) do not find any complex above the 
survey completeness limit (~ 10 5 M ) beyond a galactocentric 
radius of 4 kpc. On the other hand, star formation drops only 
at 7 kpc, well beyond the region where giant complexes are 
confined. Recent s ingle dish M33 surveys dGardan et al. 2007; 
iGratier et al.1 120 10) have shown that a population of low-mass 
molecular clouds indeed exists and becomes the dominant one 
beyond 4 kpc. Here, molecular clouds no longer aggregate 
into large complexes but form predominantly in smaller mass 
units. This is likely due to the lack of spiral arms which in this 
flocculent spiral fade away around 4 kpc. This population of 
molecular clouds may be more easily affected or dispersed by 
the growth of HII regions and therefore their properties and 
detectability might be strongly linked to the evolution of the 
associated HII region. However, since the aim of the all-disk 
surveys was to map large areas of the M33 disk, their sensitivity 
was not sufficient to unveil the presence of molecular gas around 
most of the IR sources detected by Spitzer. Thus, the nature of 
IR sources without associated Ha emission needs additional 
efforts to be clarified. 

For the aim of this paper we have restricted our sample to 
a few IR sources which are isolated, mostly located beyond 
4 kpc, in the outer regions of M33. Our sample spans a va- 
riety of F(24/im)/F(Ha) flux ratios. The IR fluxes and sizes 
of the sources suggest that they may host small young clus- 
ters or OB associations, but also evolved stars. Although there 
are no GMCs associated with these sources, the IRAM-30m 
telescope observations presented in this paper have given us 
the opportunity to investigate if these sources are associated 
to small molecular clouds (with masses < 10 5 M Q ). We have 
carried out deep searches of the 12 CO J=l-0 and J=2-l rota- 
tional lines around these sources and investigated the Spitzer 
IRAC colors of our sample and of a larger sample of AGBs 
and young clusters. There are a number of papers which in- 
vestigate diagnostic IRA C color-color diagrams to identify vari- 
ous type of sources (e.g. ICohen et al.ll2007l: IVerlev et al.l l2007: 
iGruendl & Chul 120091) . Our attempt is to distinguish in such 
color-color diagrams evolved stars from star forming regions. 

The plan of the paper is the following: in Section 2 we 
present the sample, the multiwavelength data set and the CO ob- 
servations. We discuss in Section 3 the properties of detected 
molecular clouds and in Section 4 the use of IRAC color-color 
diagrams to discriminate between AGBs and HII regions. In 
Section 5 we use the cluster birthline to investigate the massive 
stellar population in our sample. The conclusions are given in 
Section 6. 

W e adopt a distance to M33 of 840 kpc dFreedman et al.1 
1991). 



2. Molecular gas around selected mid-IR sources 

Dust emission can be investigated through the mid-IR and FIR 
data of M33 obtained with the InfraRed Array Camera (IRAC) 
and Multiband Imaging Photometer (MIPS) on board S pitzer 
(IWerneret al.l T2004: Ri eke et all [20041: iFazio et all l2004t) . The 
complete set of IRAC (3.6, 4.5, 5.8, and 8.0 fim ) and MIPS 
(24, and 70 yum) images of M33 is described in Verlev etaD 



(120071) . The image at 70 u m we use is an updated version of 
that of I Verlev etail (120071) . The total number of BCD's (Basic 
Calibrated Data) used in this new mosaic is 14980. The Mopex 
software ( Mak ovoz & Khar] 120051) version 18.3.1 was used to 
combine the data. Data are from the PID 5 "M33 mapping and 
spectroscopy" by R. Gehrz (16 AORs) and the SSC (Spitzer 
Science Center) pipeline used to produce the BCDs is the ver- 
sion S 16. 1.0. The value of the sky background removed is 9.80 
(mean value of 20 median values, each estimated in squares of 
5x5 pixels, outside the galaxy). The final measured resolution 
is 16" (gaussian FWHM * pixelsize = 13.3 * 1.2), as expected. 
In this new image the artificial stripes are less pronounced. To 
investigate the continuum ultraviolet (UV) emissio n of M 33, 
we us e Galaxy Evolution Explorer (GA LEX) dat a (Martin et al.l 
120051) . in particular those distributed bv lGil de Paz et al.l d2007l) . 
A description of GALEX observations in far-UV (FUV, 1350- 
1750 A) and near-UV (NUV, 1750-2750 A) relative to M33 
and of data reductio n and calibration procedure can be found in 
lThilkeretal.1 (120051) . 

To trace ioniz ed gas, we adopt the narrow-line Ha image of 
M 33 obtained bv lGreenawaltl (1 1 9981) . The reduction process, us- 
ing standard IRAFQ procedu res to subtract the continuum emis- 
sion, is described in detail in iHoopes & Walterbosl (l2000h . 



2.1. The sample 

We select a sampl e of 17 IR sources from the catalogue of 
IVerlev etail d2007l) . whose Ha counterparts have flux intensi- 
ties varying over two orders of magnitude. The Ha flux is gen- 
erally modest: the brightest of these sources has a luminosity 
of 1.5 xlO 36 erg s _1 which corresponds to a B0 star. We cannot 
claim there is an Ha counterpart for the faintest Ha emission as- 
sociated to a source but only some diffuse emission. In addition 
we include 3 sources located at large galactocentric radii: one 
visible only in the IR (s4), two selected in the Ha map (sl9,s20). 
In Table 1 we list the coordinates, the galactocentric distance r 
and the emission at various wavelengths of these sources. Most 
of them are barely resolved in the Spitzer map at 8 yum (2 arcsec 
resolution) and point like at 24 fim (6 arcsec resolution). In the 
Ha and FUV-GALEX maps half of them are clearly visible and 
have a small radius of 3-5 arcsec. If the AB magnitude in the 
FUV is greater than 29.5 we cannot visually identify the source 
in the map. 

The flux in the IRAC bands, at 24 yum and in Ha has been 
measured using a varyin g aperture size as d escribed in the source 
extraction algorithm of Verlev et al] (120071) . For the UV photom- 
etry we use an aperture equal to the 24 //m source size. We mea- 
sure the 70 //m fluxes in apertures which are 24 arcsec in diame- 
ter (equivalent to the FWHM size of the IRAM-30 mt CO J=l-0 
beam). The photometric errors are small, but due to the occa- 
sional presence of stripes in the images we can claim detection 
only when the flux in the aperture is above 8.5 mJy at 70 jum . 

2.2. Contamination by AGB stars 

Since 24 sources can be star forming regions or evolved 
stars (see IVerlev et al.l 12009. ), we now check whether some 
of the sources in our sample have been catalogued as evolved 
variable stars. We cross checked our source list with long pe- 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. Source sample and photometry in the MIR (Spitzer), UV (GALEX) and Ha. In column 2,3 and 4 we also give the source 
right ascension, declination and the galactocentric distance r. 



ID 


RA 


DEC 


r 


3.6 yum 


4.5 [im 


5.8 fim 


8.0 iim 


24 yum 


70 fim 


HaxlO 14 


FUV 


NUV 




deg 


deg 


arcmin 


mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


erg s -1 cirr 2 


mugAB 


mag AS 


si 


23.34048 


+30.60659 


9.20 


0.85 


0.99 


1.52 


2.85 


3.99 


(8 2) 


0.23 


30.4+3.9 


29.8+1.6 


s2 


23.35623 


+30.59784 


8.11 


0.77 


1.28 


1.54 


1.84 


1.35 




0.01 






s3 


23.55369 


+30.49387 


15.17 


1.26 


2.40 


3.11 


3.56 


3.07 




0.16 


30.6+4.2 


30.7+2.5 


s4 


23.68463 


+30.28031 


35.46 


0.62 


1.12 


2.84 


7.14 


8.70 


115 




32.4+9.6 


31.8+4.2 


s5 


23.18128 


+30.80183 


27.06 


0.14 


0.06 


0.16 


0.57 


3.93 


15.5 


0.73 


27.8+1.1 


28.2+0.8 


s6 


23.39795 


+30.43258 


14.29 


0.54 


0.35 


1.81 


5.65 


7.55 


30.7 


7.27 


27.6+1.1 


27.8+0.6 


s7 


23.43383 


+30.80995 


10.92 


0.42 


0.22 


1.81 


4.59 


4.91 


21.6 


3.16 


27.9+1.2 


28.1+0.8 


s8 


23.44272 


+30.87069 


14.59 


0.23 


0.18 


0.95 


3.18 


6.94 


25.5 


5.24 


28.7+1.7 


28.8+1.0 


s9 


23.44566 


+30.45644 


13.18 


0.62 


0.29 


0.99 


2.75 


6.19 


29.8 


21.65 


26.3+0.6 


26.4+0.3 


slO 


23.49555 


+30.92542 


16.95 


0.19 


0.14 


0.71 


2.32 


12.55 


18.3 


5.44 


28.2+1.4 


28.4+0.9 


si 1 


23 51050 


+30 96968 


19.65 


0.21 


0.19 


1.16 


2.85 


5.45 


32.5 


2.36 


28.6+1.7 


28.8+1.0 


sl2 


23.66479 


+30.52090 


21.20 


0.35 


0.16 


2.59 


9.37 


21.35 


47.7 


1.09 


29.0+2.0 


29.5+1.4 


sl3 


23.70374 


+30.62468 


20.33 


0.14 


0.08 


0.18 


0.86 


3.19 


10.9 


0.43 


31.9+7.7 




sl4 


23.24606 


+30.69297 


18.12 


3.99 


2.46 


4.21 


19.5 


12.61 


44.6 


1.56 


31.4+6.2 


30.0+1.8 


sl5 


23.55395 


+30.45810 


17.34 


0.26 


0.19 


0.38 


4.85 


14.30 


21.1 


0.82 


30.6+4.3 


30.4+2.2 


sl6 


23.56527 


+30.46270 


17.69 


1.66 


1.23 


1.30 


4.11 


11.23 


20.9 


0.45 


29.7+2.7 


29.5+1.4 


sl7 


23.61698 


+30.93344 


18.30 


0.20 


0.25 


0.52 


1.85 


7.33 


12.5 


0.28 


30.6+4.3 


31.5+3.5 


sl8 


23.69661 


+30.63216 


19.49 


0.30 


0.17 


1.23 


5.31 


5.09 


31.7 


4.57 


27.9+1.2 


28.1+0.8 


sl9 


23.67558 


+30.40020 


28.08 


0.08 


0.06 


0.30 


0.95 


0.9 


(8.1) 


1.8 


29.6+2.6 


29.4+1.4 


s20 


23.74792 


+30.61717 


23.99 


0.54 


0.71 


0.64 


0.95 






4.0 


27.8+1.2 


27.7+0.6 



riod variable stars catalogued by McOuinn et al] (120071) from 
MI R observations and w ith the variable point source catalogue 
of iHartman et alj d2006l) at optical wavelengths. We set the 
searching radius equal to the size of the 24 yum emission. We 
found that the first three infrared selected sources in our list 
(s l,s2,s3) are variables a ccording to the classification scheme 
of McOui nn et alj (120071) . The angular separation between the 
variable star and the center of the 24 yum emission is less than 
1.3 arcsec, i.e. 5 pc. We also found that the last source in our 
sampl e (s20, opticall y selec ted) is in the variable point source 
list of Hartma n et al.1 (120061) . Given the conspicuous Ha emis- 
sion associated with this region, we believe that the source s20 
is an example of an evolved star close to a young cluster (at a 
distance of about 2 pc). Given its spectral characteristics, which 
will be shown later in the paper, s4 might also be an evolved 
star: a galactic AGB or an AGB in M33 at a large galact o centri c 
distance, beyond the area surveyed by iMcOuinn et al.l ([2007). 
We keep the evolved stars in our sample to check if there is any 
detectable CO emission from the surrounding region. 

2.3. The I RAMSOmt observations 

Some CO J=l-0 emission has already been detected with a 
45 arcs ec wide beam (F CRAO) around a few sources of our 
sample dHever et al.ll2004l) . Given the large beam size, we do 
not know if the detected gas is associated to or in the proximity 
of the sources. Therefore, we have searched for CO emission us- 
ing the smaller IRAM-30 mt beam from all the sources in Table 
1. The CO J=l-0 and J=2-l lines have been observed during 
August 2007 with a FWHM beam of 24 arcsec at 1 15 GHz and 
of 12 arcsec at 230 GHz. At 24 yum all the sources are smaller in 
size than the telescope beam at 230 GHz. 

We have observed the sources in position switching mode, 
using the receiver combination A100/B100 and A230/B230 and 
the VESPA backend system with 240 MHz bandwidth (320 kHz 
resolution). One source was centered on the ON position, an- 
other source in the OFF position. The OFF source is chosen in a 
region with different line of sight velocity than the ON position 



(as seen through 21 -cm maps). The spectra have been smoothed 
to 1 km s _1 and the data from both receivers has been averaged. 

In Table 2 we summarize the CO data: integrated emission 
I (in units of main beam temperature K km s _1 ), mean veloc- 
ity V, line width W ( full width at half maximum; hereafter, 
FWHM) and peak intensity P, are estimated using gaussian fits 
to the lines. The rms refers to a spectral resolution of 2.2 and 
1.1 km s for the CO J=l-0 and J=2-l, respectively. The line 
widths have been measured using correlator spectra after cor- 
recting for hanning. 

We shall use as CO J=l-0 line intensity the average value of 
that derived by fitting a gaussian to the detected emission and 
that obtained by summing the flux in each channel inside the 
signal window (determined individually for each spectrum). The 
result is given in Table 3. As uncertainties on the CO line inten- 
sities, we consider the largest value between the following ones: 
uncertainty derived from the gaussian fit, the rms of the spectra 
integrated over the signal window, the dispersion between the 
intensity derived from gaussian fit, and the intensity derived by 
the integral inside the window signal. The gain is 6.3 and 8.7 
Jy K -1 at 1 10 and 235 GHz, respectively. Using these values, we 
derive the CO fluxes of the two rotational levels and their ratios. 
In Table 3 we also quote the H2 column densities derived from 
FCRAO (beam=45 arcsec) integrated CO J=l-0 fluxes, using a 
CO-H 2 conversion factor X C o=2.8 10 2() cirr 2 K -1 . If there are 
no entries in the column corresponding to the FCRAO column 
density, it means that the source was off the region mapped by 
FCRAO. 

The bolometric luminosity for star forming regions in Table 
3 is computed as the sum of the FUV and NUV luminosities 
uncorrected for absorption, added to the total infrared luminosity 
(hereafter, TIR luminosity): 

Lboi — Lpuv + Lnuv + Ljir + 24Lh ( 1 ) 

where the Ha lumin osity term accounts for the ionizing radiation 
dLeitherer et al.l 19991) and Ljir for the UV radiation absorbed by 
grains and re-emitted in the IR. We have not considered the con- 
tinuum radiation longward of 2800 Awhich becomes important 
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Table 2. Gaussian fits to CO lines. In column 1 and 2 we show the integrated intensity of the J=l-0 and J=2-l lines respectively, in 
column 3 and 4 the mean velocity of the two lines, in column 5 and 6 the relative line widths, in column 7 and 8 the peak intensities, 
and in the last two columns the rms for the CO J=l-0 and J=2-l respectively. 



ID 


Il-O 


I2-1 


Vl-0 


v 2 _ 


i 


w H 


Wa-i 


Pl-0 


Pw 


rmsi_o 


rms2-i 




Kkms 1 


Kkms 1 


km s" 1 


km s 


-l 


km s" 1 


km s" 1 


K 


K 


K 


K 


si 


0.284±0.028 


0.665+0.037 


-156.87+0.20 


-156. 91=1 


=0.10 


2.90+0.48 


3.19+0.31 


0.063 


0.165 


0.008 


0.019 


s2 


<0.150 


<0.200 
















0.009 


0.018 


s3 


0.247±0.025 


0.467+0.047 


-155.82+0.35 


-155.31=1 


=0.64 


5.75+1.69 


11.67+2.35 


0.039 


0.036 


0.008 


0.015 


s4 


<0.100 


<0.150 
















0.007 


0.014 


s5 


0.665±0.038 


0.746+0.107 


-179.38+0.16 


-178.93d 


=0.29 


4.93+ 0.65 


3.79+0.59 


0.109 


0.181 


0.011 


0.059 


s6 


303+0 022 

W..7V/-' _!_ V./ . \.f X_ 


639+0 038 


-108 04+0 10 


-107.89=1 


=0.08 


1.61+0.29 


1.89+0.29 


0.106 


0.222 


0.010 


0.024 


s7 


1.011±0.025 


1.338+0.033 


-240.99+0.07 


-240.94=1 


=0.06 


4.99+0.27 


4.97+0.19 


0.162 


0.244 


0.007 


0.015 


s8 


483+0 028 


981+0 035 


-243 49+0 13 


-243. 39d 


=0.07 


4 03+0 50 


3 33+0 21 


0.082 


0.235 


0.007 


0.020 


s9 


1.206±0.028 


1.921+0.035 


-124.46+0.05 


-124.36=1 


=0.03 


3.93+0.19 


3.37+0.10 


0.235 


0.494 


0.009 


0.017 


slO 


646+0 023 


1.762+0.034 


-253 66+0 08 


-253. 87=1 


=0.04 


3.61+0.22 


3 87+0 13 


0.134 


0.386 


0.008 


0.017 


sll 


1 1 14+0 027 


704+0 039 


-261 80+0 10 


-261.87=1 


-0.14 


7 64+0 65 


4 t.s+0 63 

I . *J *J -!—\J .\J^J 


0.126 


0.128 


0.007 


0.017 


sl2 


1.308±0.051 


1.516+0.042 


-151.74+0.12 


-151.40=1 


=0.06 


5.54+0.29 


5.30+0.30 


0.196 


0.290 


0.007 


0.018 


sl3 


0.809±0.024 


1.001+0.039 


-180.78+0.07 


-180.06=1 


=0.07 


3.89+0.21 


3.48+0.23 


0.163 


0.249 


0.008 


0.021 


sl4 


1.326±0.026 


1.191+0.035 


-173.56+0.05 


-173.50=1 


=0.06 


4.92+0.25 


3.90+0.24 


0.211 


0.267 


0.007 


0.017 


sl5 


0.374±0.021 


0.406+0.027 


-137.83+0.11 


-137.45=1 


=0.09 


3.11+0.66 


2.19+0.40 


0.088 


0.133 


0.008 


0.017 




0.150±0.020 


0.295+0.026 


-130.72+0.25 


-131.17=1 


=0.12 


2.33+1.45 


2.36+0.49 


0.039 


0.108 


0.008 


0.017 


sl6 


<0.100 


<0.150 
















0.007 


0.017 


sl7 


1.10±0.029 


0.838+0.058 


-261.31+0.09 


-262.45=1 


=0.16 


6.49+0.40 


4.61+0.44 


0.141 


0.171 


0.008 


0.028 


sl8 


0.975±0.030 


1.259+0.036 


-183.14+0.08 


-182.97=1 


=0.06 


5.73+0.43 


4.65+0.23 


0.148 


0.252 


0.008 


0.016 


sl9 


0.666+0.067 


0.654+0.065 


-152.91+0.49 


-154.16=1 


=0.35 


4.80+1.08 


6.32+0.91 


0.059 


0.093 


0.014 


0.028 


s20 


0.344±0.050 


0.178+0.044 


-187.52+0.57 


-190.32=1 


=0.32 


4.03+1.15 


2.03+0.86 


0.046 


0.064 


0.013 


0.029 



when young clusters have luminosities lower than 10 erg s . 
The estimated TIR luminosity, Ltir, has been computed follow- 
ing l^^etai] (Hool): 

logLr/A = logL 24 + 1 .08 + 0.5 1 x log f^L (2) 

^v,24 

and the luminosity function from the 8 and 24 fim fluxes as 
vL v . The above expression is correct for HII regions whose IR 
emission peaks longward of 24 ytzmEvolved stars sl,s2,s3 have 
no emisson in the FIR and the above formula overestimates the 
TIR. For these sources we compute the luminosity according to 
iGroenewegen et ail d2007l) who have shown that the dust emis- 
sion from evolved stars closely follows that of a black body at 
temperatures of 400-600 K. In particular, using the 8 and 24 fim 
fluxes given in Table 1, we estimate the luminosity in solar lu- 
minosity units as: 

L bol (AGB) = 130 ZL(£l) L76 e xp ( 4.4 (F 24 /F S ) 0M ) (3) 
mjy r 24 

For s4 we give the luminosities and stellar masses using the for- 
mulae for HII regions but put the values in brackets since it might 
be an evolved star. Extinction corrections for Ha fluxes in HII 
regions are 0.83 Ay where Ay, the visual extinction, is given as 

A v = 3 x0.57xlog(- — — + l) (4) 

In Ta ble 3 we list extinct ion values and cluster masses. 
Following lCorbelli et al.l d2009l) . we set 0. 1 M G as the lower limit 
of the IMF and use a Salpeter slope of 2.3 down to 0.5 M , and 
1.3 between 0.5 and 0.1 M Q . When bolometric luminosities are 
small, as in the case of these clusters, the upper end of the IMF 
is not fully populated and one has to take into account stochas- 
tic effects. In the stochastic regime, uncertainties are large and a 
simple scaling law between luminosity and cluster mass does 



not apply (and it would underestimate the cluster mass). For 
Lb i < 10 40 erg s _1 the possible range of clust er masses for 
a giv en Lt, i is up to one order of magnitude ( Corbelli et al. 
2010). In Table 3 we show the median cluster mass for the given 
bolometric luminosity, the expected median Ha luminosity, L2f 
(corresponding to the observed Lboi) and the relative uncertain- 
ties. The observed values corrected for extinction, L"J? S , are also 
given. Finally, we comment when there is not a clear peak in the 
images at the location of the 24 fim source. 

In the next Section and in the rest of this Section we discuss 
possible correlations between cloud properties. We will quote in 
parenthesis the Pearson linear correlation coefficient rp for data 
relative to star forming regions, excluding AGBs. 

The ratio of line intensities h-\lh-o in Figure 1 correlates 
with the ratio of line-widths (rp = 0.74 for log W2-1/W1-0 - 
log /2-1//1-0 )■ If we exclude sl9, the source with the largest 
W2-1/W1-0, the Pearson linear correlation coefficient is 0.91 and 
the slope is 1.9±0.2. The observed line ratio shows a marginal 
decrease as the galactocentric radius increases [rp = -0.46) and 
it does not correlate with the cluster bolometric luminosity. In 
Figure 1 we can see also that the emission at 70 //m increases 
with the cluster mass, as expected, since for most of our clus- 
ters the infrared luminosity drives the value of the bolometric 
luminosity (even though we have not used the 70 //m emission 
to compute the TIR). The scatter is however large [rp = 0.46). 
The brightness of the CO lines does not correlate with the me- 
dian cluster mass (rp = 0.38), even though Figure 1 shows that 
more massive clusters are associated with the brightest CO J=2- 
1 lines. 



3. Three methods for computing molecular masses 
and sizes 

We compute the properties of the molecular gas following 3 dif- 
ferent methods. In practice, we have three unknowns: the cloud 
size, the conversion factor Xqo between the CO luminosity and 
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Table 3. Brightness of CO clouds and associated stellar clusters. In columns 2 and 3 we give the value of the CO J=l-0 and 2-1 
line intensity averaged between that derived from gaussian fits and that obtained by summing the flux in each channel. The ratio 
between the two CO line intensities is shown in column 4. In column 5 we give the H2 column densities derived from FCRAO 
observations and in column 6 and 7 the Ha luminosity (corrected for extinction) and the bolometric luminosity of the sources. In 
the last 3 columns we show the estimated mass, expected Ha luminosity (see text for details), and visual extinction of the stellar 
clusters. 



ID 


<Ii-o > 


<h-i > 


<h-i> 


loe N'~° 

lu 5 H2FCRAO 


lop T obs 




M, 




A v 




<h-o> 






Kkms" 1 


Kkms 1 




-2 

cm 


erg s _I 


erg s 1 


M 


erg s 1 


mag 


Notes 


si 


0.300±0.028 


0.658+0.037 


2.19 


< 20.4 


36.11 


38.13" 








Visible only at 8-24-70 //m 


s2 


<0.150 


<0.200 




< 20.4 


<33.97 


37.88° 








Visible only at 8-24 //m 


s3 


0.265±0.025 


0.306+0.228 


1.15 


< 20.4 


35.99 


38.17" 








Visible only at 8-24 


s4 


<0.100 


<0.150 






34.56 


(39.00)" 


(490) 


(37.1+op 
36.4+" 
37.2 + °{ 


4.2 


Visible only at 8-24-70 fim 


s5 


0.655±0.038 


0.718+0.107 


1.10 




35.88 


38.64 


258 


0.4 


Weak Ho-.non isol. 


s6 


0.290+0.022 


0.571+0.097 


1.97 


< 20.4 


37.11 


39.15 


582 


1.0 




s7 


1.060+0.069 


1.363+0.035 


1.29 


21.02+0-™ 


36.74 


38.98 


431 


37 +o:8 


0.9 




s8 


0.565+0.115 


0.947+0.047 


1.68 


21.18!™ 


37.13 


39.00 


490 


37> 


1.4 




s9 


1.229+0.033 


1.861+0.085 


1.51 


< 20.4 


37.32 


39.30 


794 


37 4+0.6 

-Q.8 


0.3 




slO 


0.642+0.023 


1.765+0.034 


2.75 


20.90!°;^ 


37.07 


39.04 


504 


37 r 8 ° 


1.2 




sll 


1.087+0.038 


0.659+0.064 


0.61 


36.73 


38.86 


370 


36.9^f 


1.2 




sl2 


1.282+0.051 


1.681+0.283 


1.31 


20 65 +016 

— Q.26 


36.83 


39.30 


794 


-r-i 4.+O.6 

D 1 ■ -o.s 


2.3 


No FUV, tail 


sl3 


0.888+0.112 


1.000+0.039 


1.13 


9 1 9 ^ +0.05 


36.58 


38.47 


215 


36.0 + ^ 


2.7 


NoFUV 


sl4 


1.422+0.135 


1.263+0.102 


0.89 


20.78+°^ 


37.69 


39.49 


948 


37.6 + ° : * 


4.1 


No FUV 


sl5 


0.372+0.020 


0.352+0.016 


0.95 


< 20.4 


36.99 


39.11 


566 


37 - 2+ i2 

37 - 2+ n 


3.0 


No FUV, weak Ho- 


sl5 


0.159+0.020 


0.315+0.028 


1.98 


< 20.4 


36.99 


39.11 


566 


3.0 


No FUV, weak Ho- 


sl6 


<0.100 


<0.150 




< 20.4 


36.41 


38.98 


431 


37.0+™ 


2.3 


Weak Ha 


sl7 


1.159+0.082 


0.852+0.058 


0.74 


21.13 +0 11 
20.60:f| 


36.29 


38.72 


313 


36.5 + ' : » 
37.1!^ 


2.5 


Visible only at 8-24-70 jim 


sl8 


1.012+0.052 


1.269+0.036 


1.25 


36.93 


39.03 


500 


1.0 




sl9 


0.651+0.067 


0.637+0.065 


0.98 


36.48 


38.38 


206 


35.8 27 


0.9 


IR weak 


s20 


0.353+0.050 


0.172+0.044 


0.49 




36.48 


38.54 


241 


0.04 


IR weak 



" Luminosity as from equation (3), since source is an AGB star. * Source of unknown nature, luminosity computed according to equation (1). c HII 
region next to an AGB star, luminosity computed according to equation (1). 



the H2 surface density, and R the intrinsic line ratio. The latter 
quantity is the ratio between the J=2-l and J=l-0 lines we would 
measure by mapping the cloud with a telescope beam smaller 
than the cloud size. The expected line ratio is the line ratio we 
expect to observe knowing the telescope beam size at the two 
line frequencies, the intrinsic line ratio and the source extent. 
The observed line ratio is the measured value of I2-1/I1-0 If we 
consider the cloud virialized and centered at the source location, 
we have 2 equations for these 3 variables: one equation is given 
by equating the observed to the expected line ratio and the other 
one by equating the virial mass to the mass inferred from the 
CO line brightness. Clearly, even in this case some assumption 
is needed to derive the cloud properties such as mass, density, 
size etc.. We now describe 3 different methods to derive cloud 
parameters and the underlying assumptions. 

3.1. Uniform intrinsic line ratio 

We assume a fixed intrinsic line ratio R between the CO J=2- 
1 and 1-0 lines for all sources and compute the expected line 
ratio under the hypothesis that the source is isolated and centered 
in the beam. We infer the source diameter D by equating the 
expected ratio to the observed one: 



R = 



/ 2 -i(l-e^(-D 2 /(1.2®i- ) 2 )) 



/i-od 



exp{-D^ I {\.2®z-i) 2 )) 



(5) 



where is the half-width half-maximum of the beam. Knowing 
the size we can infer the virial mass and by equating the virial 
mass to the mass derived using the CO line luminosity (which 



is the same for the 2-1 and 1-0 line once we correct for beam 
dilution), we infer Xqo- Other cloud properties, such as the H2 
volume or column density, then follow. The cloud virial mass for 
a density profi le which varies inversel y to the cloud radius can 
be written as (IWilson & Scovillef^ O): 



M vh 



D , W 



— =99— (, 

M Q pc km s 1 



9 D 
f = 406 ( 



W 



arcsec km s 



(6) 



where W is the FWHM of the CO line and D is the diameter or 
cloud size. The average volume density is 



n Hl _ 406 D m ^ ms _, 2x10^ ^ 



Wf 



cm- 3 D 3 arcsec l.05 x 10 57 1.33 x m Hl 



Dl. 



(7) 



where the factor 1.33 accounts for helium. The average intrinsic 
line rati o in molecular cl ouds of our Galaxy varies between 0.5 
and 0.8 (lOka et all ! 998). being higher at low galactic latitudes. 
In exte rnal galaxies intrinsic ratios are as low as 0.3 o r as high as 
2 (e.g. lSorai etal]|2001l:lc"rosthwaite & Turnerll2007l ) and can 
be ex plained in terms of var iations of the cloud physical condi- 
tions (Sakamoto et al i ri994h . Low values are expected when the 
gas has a low density or a low kinetic temperature. In Table 4 
we give the cloud parameters derived with the first method us- 
ing a uniform intrinsic line ratio R equal to 0.5. We don't give 
the cloud size when the line ratio is larger than the maximum 
possible value (observable for a point source at the beam cen- 
ter) or smaller than the intrinsic line ratio (observable for an ex- 
tended source). For R = 0.5 estimated sizes and masses come 
out large enough that CO all disk surveys of M33 should have 
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Table 4. Cloud parameters for an intrinsic line ratio R = 0.5. The cloud diameter is given in column 2 and the virial masses 
corresponding to the J=l-0 and J=2-l CO lines are given in column 3 and 4. The H2 volume and column density are shown in 
column 5 and 6 respectively and the CO-to-H2 conversion factor in the last column (average value between those relative to the two 
CO lines). 



ID 


D 


M'r" 

vir 


M^ 1 




N„ 2 






pc 


10 4 M 


10 4 M o 


_3 

cm 


10 20 


10 20 cm" 2 K" 1 


si 














s2 














s3 


77 


26. 


107. 


15 


36 


58 


s4 














s5 


81 


20. 


12. 


10 


25 


5.6 


s6 


12 


0.31 


0.42 


54 


20 


0.7 


s7 


69 


17. 


17. 


15 


32 


4.0 


s8 


41 


6.7 


4.5 


28 


35 


2.4 


s9 


53 


8.2 


6.0 


16 


26 


1.6 


slO 














sll 


155 


89. 


29. 


7 


33 


7.3 


sl2 


57 


18. 


13. 


27 


47 


3.2 


sl3 


81 


12. 


9.7 


7 


18 


2.9 


sl4 


102 


25. 


15.6 


7 


22 


3.0 


S15 


98 


9.3 


4.6 


4 


12 


4.0 


sl5 


8 


0.52 


0.54 


168 


41 


1.0 


sl6 














sl7 


114 


48. 


24. 


9 


32 


6.0 


sl8 


73 


23. 


15.4 


19 


43 


2.8 


sl9 


94 


22. 


37. 


8 


23 


9.6 


s20 
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Fig. 1. The correlation of F70 and the cluster mass is shown 
in the upper left panel. The upper right panel shows that clus- 
ters of larger mass are associated with brighter CO J=2-l lines. 
The bottom panels display the marginal correlation between the 
line intensity ratio h-\/I\-o and the galactocentric radius and 
the good correlation between h-\lh-Q and the line width ratio 
W2-1/W1-0. Triangles are for variable stars (sl,s3) and are dis- 
played only in the lower panels. 



detected most of them. Volume densities are generally smaller 
than 10 3 cirT 3 . Since we are searching for molecular gas around 



sources one may wonder whether the assumed value of the in- 
trinsic line ratio is reasonable. The intrinsic intensity ratio can 
in fact take values significantly larger than unity when the emit- 
ting regions are warm (T > 40 K), dense (n Hl > 10 3 cirT 3 ) and 
moderately opaque (t < 5). However, larger values of R would 
imply even more massive clouds and lower densities. Thus we 
conclude that it is unlikely that the intrinsic ratio R is much larger 
than 0.5 except for the two clouds whose observed line ratio is 
larger than 2. Similarly, for s20, the intrinsic ratio is expected to 
be smaller than 0.5. 

Molecular masses according to the integrated CO line bright- 
ness can be computed as follows: 



M, 



CO 



Mr, 



41 



X 



CO 



2 x 10 2() 



D, 



I 



l-exp(-Dl rcsec /(l.2®)i 



(8) 



where the line intensity / refers to either the CO J=l-0 or J=2- 
1 line. Cloud masses derived for a constant conversion factor 
X C o = 2.8 x 10 20 crrT 2 K -1 are generally smaller than virial 
masses (except for one source in sl5). Our estimated cluster 
masses (see Table 4) are too small to account for the differ- 
ence, even taking into account the incompleteness of the IMF. 
In the last column we give X C o, normalized to 10 20 cirT 2 K _1 , 
as the value needed to have the CO-luminosity mass equal to 
virial mass (we give the average between the values relative to 
the two CO lines): 



X 



M v , 



CO 



M co 



10 



W 2 l-exp(-D 2 arcs J(l.2&, 



arc sec ) 2 ) 



D 



I 



(9) 



where the values of /, 0, W refer to the CO line which has been 
used. The large cloud sizes derived with this method for R = 0.5 
imply large cloud masses and very small star formation efficien- 
cies. Gas densities are lower than expected for typical molecular 
cloud conditions (see Section 3.4 for a comparison with Milky 
Way clouds). 
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Fig. 2. Correlations relative to parameters derived with the first 
method: in the upper panels we plot the line-width of the J=l-0 
and 2-1 line versus the cloud size. The Xco factor is plotted as a 
function of galactocentric radius and of cloud mass in lower left 
and right panels respectively. 



We do not find a correlation between the cluster mass and 
the virial mass or between the ratio M^r'/Mjr and the galacto- 
centric radius (rp = -0.15). For this model the CO-to-H2 con- 
version factor, Xco, does not correlate with galactocentric radius 
(rp = 0.34) but there is a correlation between Xco and the cloud 
mass (rp - 0.92), shown in Figure 2. This last correlation im- 
plies that Xco increases for more massive clouds. Can this be 
justified? The Xco value is expected to increase as the extinc- 
tion through the cloud or the mean volume densi ty decreases 
dMalonev & Blacklll988tlGlover & Mac Lowll2010h . Gas densi- 
ties for this model decrease as clouds get more massive, which 
justify the observed trend. Finally we find a marginal correlation 
between the cloud size and the CO linewidth, D oc W 2 , as it is 
observed in resolved molecular clouds and discussed in Section 
3.4 (>>=0.61, slope 1.9+0.6 for the J=l-01ine and r F =0.57 slope 
1.9+0.7 for the J=2-l line). 

This method is heavily dependent on the assumption that 
clouds are located at the beam center, that coincides with the 
peak of the IR-24 fim emission. Even though the molecular gas 
will likely follow the dust distribution, the peak of the mid-IR 
emission depends on the intensity of radiation field which is 
heating the grains. So the molecular cloud might be offset with 
respect to the IR emission. 



3.2. Hydrostatic equilibrium 

If cloud turbulent pressure equals the ISM ambient pressure 
given by hydrostatic equilibrium we can estimate the gas vol- 
ume density of the cloud. This is in reality a lower limit since 
cloud cores can be denser and self gravitating, and the outer parts 
close to equilibrium with the surrounding ISM. From the mini- 
mum gas density we infer D max , the maximum cloud size, using 
the virial equation. As we shall remark in Section 3.4, it might be 



that only part of the cloud is gravitationally bound and in virial 
equilibrium. In this case the virial equation applies only for the 
higher density core with size D < D max . 

In Table 5 we give the resulting cloud parameters and below 
some details of this model. We compute the ambient pressure as 



Pism = ^GI, g (E g + ^£ s ) cm 3 K 



(10) 



(ICorbellil f2003) where E gs are the gas and stellar surface mass 
densities and c g)S are the velocity dispersion relative to gas and 
stellar disk. The velocity dispersion of the gaseous disk is c g = 
8 km s _1 . We shall use the local values of the HI column densities 
for the gas surface densities (WSRT data) and a scale height of 
the stellar disk zo = 0.5 kpc for the stellar dispersion. The stellar 
dispersion and surface densit y derived thro ugh the dynamical 
analysis of the rotation curve (ICorbellill2003l) reads: 



c s = (tiGi^s) 



S s = 430 exp(-R(kpc)l\A2) M e /pc 2 (11) 



The ISM pressure in M33 does not vary much radially since 
the gas surface density has a very large scale length. However it 
can vary locally due to the presence of filaments and flocculent 
spiral arms. For each source the turbulent pressure is 



1.16n H ,r ; 



Turb 



25 n H , W 2 (km/s) cm' 3 K 



(12) 



(1.16 accounts for Helium, see also Spitzer 1978, p. 37). The gas 
volume density is then: 



/Wcm- 3 K) 
25 W 2 (km/s) 



(13) 



We shall use the smallest linewidth measured in the 1-0 or 
2-1 line, which gives the highest volume density. We compute 
the maximum cloud size by equating the virial mass to the mass 
in pressure equilibrium which reads: 



l74W 2 (km/s) W 2 (km/s) 

D arc sec = *l ~ = 66 " 



n 



//, 



(14) 



ISM 



This relation naturally explains the size-linewidth relation 
W oc ^^D. In fact we find that the 2-1 linewidth (which is mostly 
used in this relation) correlates with the source size, with the 
right slope (see Figure 3, rp = 0.94 slope=2.0+0.2). A similar 
relation exists for the 1-0, but with a larger scatter (rp = 0.81). 
The size-linewidth relation implies another important feature of- 
ten mentioned about molecular clouds, namely their rather con- 
stant column densities. For this model in fact the gas column 
density reads: 



7V™"(cm- 



10 20 



0.1 26D" 



c n™(cm- 3 ) oc Vi 7 



ISM 



(15) 



Values of the local ISM pressure in M33 do not vary much 
and this implies almost constant column densities, as can be seen 
in Table 5. We can also infer the maximum value of Xco from the 
virial mass to the CO-luminosity mass ratio using the maximum 
size D: 



vmax 

A CO _ 



My 

Men 



= 10 



W 2 l-exp(-(D'^ ec ) 2 /(l.2®) 

jynxax J 



(16) 



Where the values of /, 0, W refer to the line which has 
been used (smaller linewidth). The correlation between cloud 
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Table 5. Clouds in pressure equilibrium with the surrounding ISM. The CO rotational line used to infer the cloud parameters is 
given in column 2, the maximum value of the intrinsic CO line ratio in column 3 and the maximum allowed cloud diameter in 
column 4. The maximum cloud mass and the corresponding minimum volume and column density of H2 are given in column 5,6 
and 7 respectively. The maximum value of the CO-to-H2 conversion factor is shown in the last column. 



ID 


line 




pc 


M max 
10 4 M o 


Hi 

cirT 3 


fjmm 
Hi 

10 20 


vmax 

10 20 


si 


1-0 


0.66 


41 


3.5 


14. 


18 


6.4 


s2 
















«3 


1 n 


u.uu 




j j . 


8 fS 


78 




s4 


















2- 1 




04 


14 


4 fS 




8 9 
0. j 




1 n 


SO 


1 
1 u 




80 


24 


n s 




2-1 


U.JO 


81 




1 ] 


78 


8 9 
0. j 




2-1 




41 


4 7 


1 8 




4 fS 


s9 


2-1 


0.43 


36 


4. 1 


26 


29 


1.9 




1 n 

i-u 


u.oo 


4Q 


f\ 4 


1 5 

1 j 


Z. J 


S 4 


sll 


2-1 


0.22 


61 


5.9 


15 


28 


5.5 


sl2 


2-1 


0.63 


90 


25. 


11 


30 


8.6 


sl3 


2-1 


0.39 


57 


6.9 


11 


19 


4.0 


sl4 


2-1 


0.27 


41 


6.3 


25 


31 


2.5 


sl5 


2-1 


0.24 


13 


0.6 


85 


33 


1.0 


sl5 


1-0 


0.51 


14 


0.8 


76 


33 


2.9 


sl6 
















sl7 


2-1 


0.31 


77 


16. 


11 


26 


6.7 


sl8 


2-1 


0.47 


65 


14. 


15 


30 


7.0 


sl9 


1-0 


0.77 


147 


34. 


3.1 


14 


15. 


s20 


2-1 


0.13 


23 


0.9 


23. 


16 


1.5 



mass and the CO-to-H2 conversion factor is shown in Figure 3 
(rp = 0.95). From the maximum size we can infer the maximum 
intrinsic line ratio, R max , given in Table 5 for isolated sources, at 
the beam center. 

With this method cloud sizes and masses are less dependent 
on the assumption of the cloud being centered with respect to 
the millimeter telescope beam. Unfortunately we are able only 
to derive upper limits to M and D. In Section 3.4 we use the ob- 
served size-linewidth relation to infer the size of self-gravitating 
cloud cores. We notice that there are 3 sources with X" Q X < 1.5. 
Those are sources with sizes smaller than 25 pc and small masses 
(< 10 4 M ). Since M33 has a lower metal abundance than the 
Milky Way we do not expect an overabundance of CO. Hence, 
these sources are proba bly overlum inous in CO or they are not 
in virial equilibrium (Malonevl l 19901) . The molecular gas might 
be so warm that the J=l-0 luminosity is no longer proportional 
to the H2 column density. For example for a gas at 40 K the J=l- 
lin e intensity is more th an twice as bright as a gas at 20 K 
(e.g. lSakamoto et aUll994l) . Recent high resolution observations 
in the M3 3 outer disk have d iscovered clouds with high CO lu- 
minosity (Bigi el et al.ll20 10) in the proximity of sources host- 
ing massive stars. The very small values of R which we recover 
for some sources might be indicative of an offcenter molecular 
cloud due to the presence of massive stars. This likely holds for 
sl5 where we detect two components, and for s20 where the CO 
observations have been centered on the optical HII region in the 
absence of a 24 /vm counterpart. 

3.3. Radially varying X co 

In Table 6 we derive the cloud properties for a radially varying 
Xqo factor. The radial scaling used takes into account t he we ak 
Xqo metallicity dependence estimated bv lBolatto et al.l (12008): 

logXco =a- 0.35 x log- => X co = 2.8 x io 20+003r (17) 



where the average metallicity gradient O/H used is given by 
Magrini et al. (2010) and the galactocentric distance r is in kpc. 
The value of the constant a is such that Xco is 2.8 10 20 cm -2 K _1 
at the center of M33. The metallicity dependence of the Xco fac- 
tor has not yet been firmly established. We prefer to consider a 
shallow radial dependence because in M33 the radial decrease of 
the average radiation field, which dissociates the CO molecule, 
might balance the already shallow radial metallicity gradient and 
keep Xco almost constant. 

We derive the source size by equating the virial to the molec- 
ular mass given by the CO line luminosity: 

M * 4Q6 Wl l-eM-DL s J(l-2®) 2 
M co ^IXco I D 

From the size we infer the volume density. We then use the 
source size and the observed CO line ratio to infer the intrinsic 
line ratio. We take the CO line with the largest linewidth because 
it gives the highest gas density (see Eq.(6)).EI 

For isolated clouds, centered in the beam, the intrinsic line 
ratio R expected for that size is also shown in Figure 3 and de- 
creases moving radially outwards. 

Limits to the validity of this method come from the fact that 
the weak metallicity dependence used here for Xco is based 
on CO interferometric measurements which might underesti- 
mate cloud sizes a nd hence invalidate the scaling law used 
(llsrael & B aas 2001). Metallicity in M33 however, has a scatter 
at each radius larger than the overall observed radial gradient: 
a stronger metallicity dependence can be considered only if the 
metallicity at each source location is known. This method has 
been based on the assumption that clouds are in virial equilib- 
rium. In our galaxy data suggest that clouds with masses below 

2 We call some caution for this choice due to the possibility of de- 
tecting more than one cloud in the source proximity, especially since it 
is often the J= 1-0 line, observed with the largest beam, which has the 
largest linewidth. 
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Table 6. Cloud parameters for a radially varying Xqo- The CO rotational line used to infer the cloud parameters is given in column 
2, the intrinsic CO line ratio in column 3 and the cloud diameter in column 4. The cloud mass, the H2 volume and column density 
are shown in column 5, 6 and 7 respectively. 



ID 


line 


R 


D 


M 


l-U 
Hi 


Nj/ 2 








pc 


10 4 M Q 


cm 4 


10 20 cnr 2 


si 


2-1 


0.566 


17 


1.7 


98 


51 


s2 














s3 


2-1 


0.289 


1.2 


1.6 


270000 


9990 


s4 














s5 


1-0 


0.287 


21 


5.0 


160 


104 


s6 


2-1 


0.670 


56 


2.0 


3 


5.2 


s7 


1-0 


0.344 


26 


6.3 


112 


90 


s8 


1-0 


0.441 


33 


3.5 


98 


100 


s9 


1-0 


0.501 


53 


8.1 


16 


26 


slO 


2-1 


0.749 


29 


4.2 


53 


47 


sll 


1-0 


0.154 


13 


7.3 


1058 


424 


sl2 


1-0 


0.360 


33 


9.1 


100 


102 


sl3 


1-0 


0.340 


43 


6.4 


24 


32 


sl4 


1-0 


0.264 


41 


9.9 


42 


53 


sl5 


1-0 


0.254 


26 


2.5 


43 


35 


sl5 


2-1 


0.515 


19 


1.0 


46 


27 


sl6 














sl7 


1-0 


0.191 


19 


7.7 


360 


211 


sl8 


1-0 


0.329 


21 


6.9 


213 


138 


sl9 


2-1 


0.249 


13 


5.1 


713 


286 


s20 


1-0 


0.125 


16 


2.6 


185 


91 
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Fig. 3. Correlations relative to parameters derived with the sec- 
ond method (a,b,c) and with the third method (d) 



3.4. Comparison with spatially resolved molecular clouds 

Recen t papers (e.g. iBolatto et al.l l2008t iRoman-Duval et aTl 
2010) summarize relevant physical properties of Milky Way and 
extragalactic molecular clouds. Masses of clouds inside the solar 
radius mostly range between 10 2 and 5 x 10 5 M and their sizes 
are below 25 pc. The median valu e of the average H2 volum e 
density distribution (see Table 1 of IRoman-Duval et aT]|2010l ) 
is 230 cirT 3 with half of the clouds having densities between 100 
and 400 cm -3 and with low mass clouds having the highest den- 
sities. The BIMA all-disk survey of M 33 has identified GMCs 
with a completeness lim it of 10 5 M (lEngargiola et al.l 2003: 
iRosolowskv et al.ll2.007l) . Hence, we expect our sources to be as- 
sociated to smaller mass clouds, below the completeness limit of 
the BIMA survey. M33 is a low-luminosity galaxy with a lower 
gas and stellar surface density than the Milky Way and flocculent 
spiral arms. The mass spectrum of molecular clouds in M33 is 
steeper than that of the Milky Way dBlitz & RosolowskvlkOOSl) 
and the average cloud properties are more similar to those of the 
clouds in the outer Galaxy. Molecular clouds in fact respond to 
large scale environmental variations. In the outer Galaxy molec- 
ular clouds have smaller masses than in the inner regions and 
it is still a debated question whether clouds of small mass are 
pressure bounded or in virial equilibrium dHever et aTl )200 1 ; 
Blit z & Rosolowskvll2005tlHever et alj|2009b . 

For a quantitative comparison between cloud properties in 
our sample, inferred using the three models presented in this 
Section, and resolved molecular clouds, we shall focus on the 
observed size-linewidth relation: 



10 4 M are strongly influenced by nongravitati onal forces an d 
virial equilibrium overestimates their masses (Maloney 1990). 
Using only the CO observations presented in this paper it is hard 
to establish if the virial equilibrium assumption is correct for 
clouds in our sample. If not, clouds might be pressure confined 
(as considered in Section 3.2) or related to short-lived event (e.g. 
interstellar shocks). 



W= 1.03 XD' 



0.5 



(19) 



where W is the FWHM of the CO line in km s and D is the 
cloud diameter in p arsec. The exponent of the relation was orig- 
inally derived by ISolomon et al.1 (119871) . while the coefficient 
1.03 comes from the best fit to resolved clouds in the Milky 
Way and in M33 (IRosolowskv et al.ll2003t iBlitz & Rosolowskvl 



10 



Edvige Corbelli et al.: On the nature of faint mid-infrared sources in M33 



I2005h . We use 30% larger cloud radii for M33 data to over- 
come the difference in the size m easurement technique wi th the 
Milky Way data pointed out by iRosolowskv et a l. (2003). The 
above relation fits the data over two orders of magnitude in ra- 
dius, from 1 to 100 pc, and it is shown by the continuous line 
in the three panels of Figure [4] In the same Figure we plot the 
data for M33 GMCs (op en symbols) dWilson & Scoviilelll990t 
IRosolowskv et al.l 12003) and for clouds relative to HII regions 
in our survey (filled circles). Cloud sizes derived using a uni- 
form intrinsic line ratio R=0.5 (bottom panel) are too large to 
satisfy the observed size-linewidth relation. Furthermore, cloud 
masses are too large to be compatible with previous surveys and 
the median volume density is lower by more than a factor 10 
than the value in the Milky Way. Considering the outer envelope 
of the clouds in hydrostatic equilibrium with the ISM, we esti- 
mate higher densities. However we cannot derive the parameters 
of the cloud using only this assumption but their limiting values. 
To overcome this limitation, we assume that the size-linewidth 
relation observed in the Milky Way and in the GMCs of M33 
holds for gravitationally bound cloud cores, and that there is an 
average scaling factor of the core size with respect to the max- 
imum cloud size jy max (given in Table 5 and plotted with aster- 
isk symbols in the middle panel). By minimizing the dispersion 
around the plotted size-linewidth relation, we derive core sizes 
which are a factor 4 smaller than D" ,£ ". These are shown in the 
middle panel of Figure|4]by filled circles. In the top panel of the 
figure we plot cloud sizes relative to the radially varying Xco 
model (Table 6). We can see that despite the large scatter the 
average cloud size is smaller than that of GMCs and not much 
different than what the size-linewidth relation predicts. 

4. IRAC colors and the search for embedded 
sources 

In the previous Section we have seen that we were able to de- 
tect CO in most star forming regions and in the proximity of 
some evolved stars, even though CO lines are generally weaker 
in the latter case. We shall now examine whether newly born 
clusters can be distinguished from evolved stars using IRAC 
color-color diagrams. Previous studies on closer objects (like 
the LMC and SMC galaxies) used IRAC color-color diagrams 
to separate Young Stellar Objects (hereafter YSOs) from stars or 
from evolved AGBs, and AGBs from HII regio ns (Simon et al. 
20071: ICohen et al.l l2007t iMarengo et al] 120081: [Buchanan etai] 



2009; lGruendl& Chul 120091) . Models indicate that star forma- 
tion regions at different stages of the protostar to star evolution- 
ary sequence occupy different areas of the color-color diagrams, 
even though there is overlappi ng and not a well defin ed separa- 
tion between the variou s areas dRobitaille et al.1120 06). Similarly 
iGruendl & Chul (120091) have shown that there is no simple di- 
agnostic in color-color or color-magnitude diagrams that can be 
used to uniquely separate YSOs from AGB stars. Evolved stars 
of different types (Carbon-rich versus Oxygen-rich for exam- 
ple) lie in separate regions of the IRAC color-color diagram bu t 
HII regions are found to lie everywhere (Bucha nan et alJ l2009). 
For AGB stars the [3.6]-[8.0] and [3.5]-[4.5] colors are used to 
separate carbon stars from stars with silicate envelopes due to 
the presence of silicate and H2O features in the spectra. These 
feature affect the [3.6]-[8.0] color excess which seems to corre- 
late with the mass-loss rate of individual stars and the optical 
thickn ess of their envelope (Marengo et al. 2008: iGroenew egen 

We display the mid-IR SEDs between 3.6 and 24 /mi in 
Figures 5 and 6. In Figure 5 we select all sources where AFx 




og D/2 (pc) 



Fig. 4. The solid line is the fit to the size-linewidth rela- 
tion for the Milky Way and M33 spatially resolved molec- 
ular clou ds. The open sym b ols ar e for GMCs observed in 
M33 by Rosolowskv et al] (120031) (open squares) and by 
IWilson & Scovilld (11990) (open triangles). The filled circles are 
molecular clouds associated to our HII regions. Sizes have been 
computed according to the uniform intrinsic line ratio in the bot- 
tom panel and to the radially varying Xco model in the top panel. 
In the middle panel asterisks indicate the maximum cloud size 
for hydrostatic equilibrium models and filled circles the core 
sizes, which are 0.25 the maximum size. This scaling minimizes 
the scatter around the size-linewidth relation. The linewidths 
used are the CO J=l-0 linewidth in the bottom panel and are 
given in Table 5 and Table 6 for the other two panels. 



varies by less than one order of magnitude between 5.8 and 8 /mi 
, while sources with larger differences are plotted in Figure 6. 
Notic e that all sources in Figure [6] have a dip at 4.5 /mi (see 
also IVerlev et al.1 120071) except si 7, which is the source with 
the lowest Ha flux. Sources are ordered according to the SED 
slope between 8 and 24 /mi with s5 being the source with 
the largest gradient and sl4 the source with the smallest one. 
For all sources in Figure [6] we detect CO emission. Sources in 
Figure [5] instead have either weak or no CO emission and no 
dip at 4.5 /mi. We have identified sl,s2,s3, and s20 as regions 
containing evolved variable stars. It has been shown that the 
MIR SED of YSOs and HII region s have steeper slopes t han 
evolved stars dBuchanan et al.l l2009t IGruendl & Chull2009l) . in 
agreement with our finding. This is because the SED of star 
forming regions i n the mid-infrared pe aks at wavelenghts longer 
than 24 /mi (e.g. lLawton et al.l 2010, ). The CO emission as- 
sociated with sources in Figure [5] is weak or under our detec- 
tion limit. The McOuinn et al.l d2007l) catalogue does not extend 
as far as s4, which might be a variable associated with M33 or 
with our Galaxy. The only HII region where no CO has been 
detected is si 6, likely an evolved HII region. Evolved stars and 
evolved HII regions show a shallow slope in the IRAC bands. 
The dip at 4.5 /mi is present in the SEDs of HII regions even 
though some star forming site (si 7 for example), might not have 
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Fig. 5. IRAC SEDs: sources with a small flux increase between 
4.5 and at 8 //m. The plot shows the sources from bottom to top 
scaling the fluxes arbitrarily so that the source at the bottom is 
the one with the largest gradient between the 8 and the 24 fim 
flux and the source at the top is the one with the smallest gra- 
dient. Effective flux values can be recovered from Table 1 . All 
sources except sl6 show a monotonic increase of the mid-IR 
flux between 3.6 and 24 /mi. The 24 fim flux for s20 is below 
threshold. 



in 



o 




Fig. 6. IRAC SEDs: sources with a large flux growth between 
4.5 and at 8 fim. The plot shows the sources from bottom to 
top scaling the fluxes arbitrarily so that the source at the bottom 
is the one with the largest slope between the 8 and the 24 /mi 
flux, the source at the top is the one with the smallest gradient. 
Effective flux values can be recovered from Table 1 . Open dots 
are used for sources which do not show a noticeable decrease of 
the flux between 3.6 and 4.5 /mi. 



it. Some YSO SED s show that the dip is common b ut not al- 
ways present dGruendl & Chdl2009t Ishnon et al.ll2007l) . Results 
from ISO spectroscopy around massive protostars have shown 
that C O and CO? absor ption bands are present in the 4.5 /mi 
band ( Ivan Dishoeckll2004r) . These features can explain the dip 
mostly present in the SED of star forming regions. The CO2 ab- 
sorption line and a time dependent CO ab sorption feature can 
also be present in the s pectra of AGB stars (Marengo et al.l2010t 
Justtanont et al. 2004) but their strengths are not sufficient to af- 
fect notably the flux in our sample. 

Fig. [7] shows al l long period variables catalogued by 
iMcOuinn et al.l d2007h which have been detected at 3.6, 4.5 and 
8.0 fim in the [3.6]-[8.0]/[3.6]-[4.5] and [4.5]-[8.0]/[3.6]-[4.5] 
color-color d iagrams. In the sam e figure we plot sources in the 
catalogue of Verle v et all (120071) with Ha emission, and hence 
are likely to be young star forming regions, and do not have any 
long-term variable coincident with the extent of the 24 fim emis- 
sion. As we can see, variable stars occupy a well defined area and 
HII regions have a distribution which is centered on a different 
region of the color-color diagram, with some overlap. 

The two lines in Figure|7]have slopes 1 and 0.5 respectively 
and separate the region populated by evolved stars. Only about 
3% of the evolved stars lie to the right of the lines (about 10% 
of the IR selected HII regions lie to their left). The sources con- 
sidered in this paper well represent the two distributions. Our 
HII regions with a detectable CO, 24 fim and Ha emission lie 
in a well defined region and not in the upper left corner where 
most of the evolved stars lie. Close to the dividing line we find 
si 7, an embedded HII region candidate without Ha and very low 



UV fluxes, and si 6, the HII region devoid of CO. We cannot ex- 
clude that si 6, sl7 host an evolved star which did not meet the 
classification schemes used (actually they are bo th within a few 
parsec of a non-point source variable according to Hartma n et al.l 
(l2006h ). The probability for an AGB s tar to encounter a SF re- 
gion is low in the solar neighborhood (Kastner & Mvers 1994) 
but it can be appreciable along spiral arms and filaments. The 
colours of s20, an evolved HII region with no 24 fim counter- 
part and associated to an evolved variable star has IRAC colors 
similar to variable stars. 

In Figure [8] we display other variables and colors, even out- 
side the IRAC range. In each panel the young HII regions where 
dust and CO emission have been detected occupy well defined 
areas. For example all of them except one (sl4) have the 3.6 /mi 
magnitude above -1 and have [24]-[3.6]< -6.5. Similarly their 
IRAC color [5.8]-[8.0] is > 1.6. There is also a nice correlation 
between the [4.5]-[5.8] color and the [3.6]-[8.0] color with the 
young HII regions where dust and CO emission have been de- 
tected lying at [4.5]-[5.8]> 1. Triangles in the upper panels of 
Figure [9] i ndica t e the colors of 24 fim sources in the catalogue 
of I Verlev et al.l (120071) with no Ha counterpart which "host" a 
variable star i.e. at least one classified variable is located within 
the 24 fim e mission boundary. O nly 25 out of 125 sources in the 
catalogue of I Verlev et al.l d2007l) with no Ha counterpart and de- 
tected in IRAC bands have been classified as variable stars. Most 
of these lie to the left of the lines drawn in Figure [9] Variable 
sources which lie to the right of the line may contain both an 
evolved star and an embedded HII region. In the bottom panels 
we show sources with no Ha counterpart according to the cat- 
alogue of IVerlev et al.l (120071) not associated to known variable 
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Fig. 7. Bottom panels display in the IRAC color-colo r diagrams 
the variable stars detected by McQuinn et al] d2007l) (small dot 
symbols) and young HII regions (the open circles; these are 
sources in the catalogue of lVerlev et al.l d2007f) with Ha counter- 
part and not classified as variables) . Only 3% of variable stars lie 
to the right of the continuous line. In the upper panels we show 
the same IRAC colors but for our source sample. Open symbols 
are sources where no CO has been detected, filled symbols are 
sources with CO detection. Triangles are for variable stars. The 
open circle is for source 4, with uncertain classification. Squares 
are for sources from 5 to 19, i.e. HII regions with 24 //m coun- 
terparts. The filled dot is for s20, a variable star coincident with 
an HII region without 24 /mi emission. 
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Fig. 8. Other colours... same symbols as in previous Figure. The 
source with the highest flux at 3.6 fim is sl4. 
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Fig. 9. Star symbols in the lower panels are for 24 yum sources 
with no Ha co unterp art which have been not identified by 
iMcOuinn et al.l (120071) as variable stars . Filled triangles in the 
upper panels are sources in Verle v et al.l d2007l) catalo gue which 
host variable stars according to McQui nn et al.l d2007l) . 



stars. The number of these sources which lie to the left of the 
line exceed by far the percentage expected if these sources have 
similar colors to the visible HII regions. Either embedded stars 
occupy a wider area of the IRAC color-color diagram than HII 
regions or there are unidentified AGBs, variables in crowded re - 
gions which escaped the classification by McQui nn et alJd2007l) . 
The majority of sources to the right of the line are good candidate 
for being young star forming regions. Only 5 of these sources are 
associated with known GMCs. The associated molecular clouds 
might be of smaller mass, like the ones detected in the CO survey 
presented here. 

5. The cluster birthline test 

To verify the young age of our sources we will check if they 
lie close to the cluster birthline, which is the lin e in the Li, n i- 
Lboi/L Ha plane around which young clusters lie dCorbelli et al.l 
2009). Older clusters should be located above the birthline be- 
cause of the faster fading of the Ha luminosity compared to the 
bolometric luminosity. Leakage of ionizing photons or the ejec- 
tion or the delayed formation of massive stars can also move 
the clusters above the birthline. In Figure [10] we show with star 
symbols the birthline relative to a stochastically sampled uni- 
versal IMF with a Salpeter slope -2.35 at the high-mass end. 
Stochasticity allows a spread in the values of ~Lboil^Ha for a given 
Lt,oi: for a given bolometric luminosity we can have a fully pop- 
ulated IMF up to a certain mass M max , or an IMF populated only 
up to M, < M max , plus a brighter outlier. The errorbars in each 
bin indicate the dispersion due to bright outliers that can form for 
the given L^ /. The highest simulated L/, //Lhq. value is for a fully 
populated IMF case. We see that there are two sources (sl2,sl6) 
which are clearly above the upper 3-cr boundary. These might 
be more evolved sources or leaking ionizing radiation. Since we 
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Fig. 10. Sources on the cluster birthline. Star symbols indicate 
the birthline relative to a stochastically sampled universal IMF 
with their 1-cr dispersion. Open symbols are for sources where 
no CO has been detected, filled symbols are sources with a CO 
detection. The filled dot is for s20 which does not have IR emis- 
sion. 

have not found CO emission from si 6, it is likely that it is an 
evolving source. The source sl2 is the brightest in the IR and 
has a strong CO emission: thus, its position can be explained if 
it is in the process of forming its massive stellar population or if 
it is leaking ionizing photons. The strong 8 /mi emission places 
sl2 in the upper right corner of the [3.6]-[8.0]-[4.5]-[5.8] color 
plot. All the other sources cannot be found below the birthline if 
the IMF assumptions are correct. Our sources are all compatible 
with a random sampled IMF model including the two sources, 
sl3 and si 9, close to the lower boundary which might host out- 
lier stars. 

6. Conclusions 

We have investigated the nature of faint 24 //m sources with a 
weak or no Ha counterpart. These could be regions with embed- 
ded clusters still associated to molecular gas, or evolved stars 
with dusty envelopes. To distinguish between these two possi- 
bilities, we have used published catalogues of variable stars in 
M33 and carried out deep observations of the CO J=l-0 and 2-1 
line with the IRAM 30-mt telescope to search for molecular gas 
associated with 20 sources. The main results are: 

- Deep pointed CO observations have revealed the presence 
of CO around 17 of the 20 sources included in our sam- 
ple. The weakness of the CO lines around most of our 
sources is indicative of the existence of a large population 
of faint CO clouds in M33. This is in ag r eemen t with the 
results of the CO survey of Gra tier et alj ([2010) and with 
the steeper molecular ma ss spectrum relative to our Galaxy 
(Blitz & Rosolowskvll2005l) . 

- The estimate of cloud masses is uncertainsince our pointed 
observations are deeper than any existing CO map. 



Considering clouds in virial equilibrium and a weak ra- 
dial dependence of the CO-to-H2 conversion factor, we find 
molecular gas masses in the range 10 4 - 10 5 M . In this 
case, cloud sizes agree on average with those predicted by 
the size-linewidth relation of the Milky Way and resolved 
GMCs in other galaxies, but have large dispersions. We have 
also estimated limits to molecular gas masses and sizes in 
the case where cloud outer envelopes are in pressure equilib- 
rium with the surrounding medium. This model reproduces 
the observed size-linewidth relation correctly if gravitation- 
ally bound cloud cores are a factor 4 smaller in size than the 
outer envelopes. 

- The intrinsic CO J=2-l to J=l-0 line ratio is generally small, 
~0.4, if one makes the assumption that clouds are within a 
few parsecs from the beam center. The CO lines are weaker 
at the location of evolved stars than around star forming re- 
gions, as is the Ha and 70 /mi emission. Stellar clusters asso- 
ciated to clouds in our sample are young since they lie along 
the cluster birthline. 

- In the IRAC color-color diagrams AGB variable stars and 
HII regions occupy distinct areas even though there is some 
overlap: most of the AGBs have [3.6]-[4.5] color in the range 
to 1.5 and lie around a line of slope 0.5 in the [3.6]-[8.0]- 
[3.6]-[4.5] plane. Most of the 24 /mi sources associated to 
HII regions have higher values for the [4.5]-[8.0] or [3.6]- 
[8.0] color than variable stars. The IRAC colors of our 20 se- 
lected sources are consistent with these distributions. They 
show a linear correlation between the [4.5]-[5.8] color and 
the [3.6]-[8.0] color with the AGBs having the lowest val- 
ues. The IRAC colors confirm the presence of AGBs in 4 
sources of our sample whose variability was previously de- 
tected either in infrared or optical surveys. 

- Using IRAC color-color diagrams we also predict that in the 
sample of 24 /mi sources without Ha counterpart there are 
still unidentified evolved stars other than young clusters. The 
young clusters might lack Ha emission because they are of 
small mass or because they are very compact, and deeply em- 
bedded. Additional sensitive CO observations and high res- 
olution far-infrared images will help to identify unambigu- 
ously the properties of the young star forming sites in M33. 
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